arXiv:1502.06059v2 [cond-mat.mtrl-sci] 2 Apr 2015 


H-Si bonding-induced unusual electronic properties of 
silicene: a method to identify hydrogen concentration 


Shih-Yang Lin^, Shen-Lin Chang^’^’*, Ngoc Thanh Thuy Tran^’^, Po-Hua Yang^, and Ming-Fa Lin^’^ 
^Department of Physics, National Cheng-Kung University, Tainan 701, Taiwan 
^Department of Physics, National Sun Yat-Sen University, Taiwan 
^National Center for High-Performance Computing, Taiwan, 710 Tainan, Taiwan 

April 3, 2015 

Abstract 

Hydrogenated silicenes possess peculiar properties owing to the strong H-Si bonds, 
as revealed by an investigation using first principles calculations. The various charge 
distributions, bond lengths, energy bands, and densities of states strongly depend 
on different hydrogen configurations and concentrations. The competition of strong 
H-Si bondings and weak sp^ hybridization dominate the electronic properties. Chair 
configurations belong to semiconductors, while the top configurations show a nearly 
dispersionless energy band at the Fermi level. Both two systems display H-related 
partially flat bands at middle energy, and recovery of low-lying vr bands during the re¬ 
duction of concentration. Their densities of states exhibit prominent peaks at middle 
energy, and the top systems have a delta-funtion-like peak at i? = 0. The intensity 
of these peaks are gradually weakened as the concentration decreases, providing an 
effective method to identify the H-concentration in scanning tunneling spectroscopy 
experiments. 
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1. Introduction 


Two-dimensional materials composed by group IV elements have attracted consider¬ 
able attention in the helds of chemistry, material science and physics,^”® mainly owing to 
their planar structure with sp^ or sp^ bondings formed by four valence electrons. This 
has provided the motivation to conduct studies on their electronic,optical,and 
transportproperties. The gap opening in the electronic properties is an important 
issue for its possible use in future nanoelectronic devices.Several experimental meth¬ 
ods have been realized, one of which, the chemical functionalization,^®”^^ is particularly 
promising. Recently, the successful hydrogenation of graphene has created a novel material 
known as graphane with a large band gap,®^”^® which has been verified by angle-resolved 
photoemission spectroscopy (ARPES).®^”^® During the synthesis process, its band struc¬ 
ture can be signihcantly altered and controlled through reversible hydrogenation,®^ while 
the concentration is very difficult to estimate. However, graphene faces many challenges 
in applications, such as toxicity, difficulty in processing, issues to integrate it with current 
silicon-based electronic technology, and the availability of large-area sheets. This work in¬ 
vestigates in detail how the unique H-Si bonds with special charge distributions induce the 
feature-rich geometric and electronic properties in hydrogenated silicene. 

As the counterpart of graphene, 2D hexagonal lattices of Si, namely silicene, and Si 
nanoribbons have been fabricated by the chemical vapor deposition process on a silver 
substrate®’®®”®® and Iridium substrate.®® Other chemical exfoliation methods have suc¬ 
ceeded in creating suspended silicon monolayer sheets doped with magnesium and covered 
with organic groups.These experimental realizations of providing suitable environ¬ 
ments for single-sided and double-sided chemical functionalizations further tune the elec- 
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tronic properties of silicene drastically. Quite different from graphene or other carbon-based 
nanostructures, silicene with its diverse electronic properties could enable the electronic in¬ 
dustry to produce fast nanoscale electronics without the need to retool. 

Unlike graphene, silicene is not stable as a perfectly planar sheet. Because the sp^ hy¬ 
bridization is more stable in silicene than the sp^ hybridization, it is energetically favorable 
as a low buckled structure, which has been conhrmed by DFT studies and phonon disper¬ 
sion calculations."^^ Recently, the electronic structure of silicene has been revealed to be 
that of a zero-gap semiconductor by ARPES and scanning tunneling spectroscopy (STS) 
measurements.Other theoretical works exploring silicene sheets passivated with hydro¬ 
gen show the emergence of an indirect middle band gap.^®’"^® These results suggest that the 
electronic properties of hydrogenated silicene can be altered by a competition between sp^ 
hybridization and hydrogenation. However, since silicene sheets have only been recently 
realized, the effects of various hydrogen concentrations with different conhgurations on this 
new material have yet to be thoroughly explored. 

In this paper, how to create the diversihed electronic properties via competition between 
the H-Si bonds and the Si-Si a and vr bonds is investigated by hrst principle calculations. 
The double-sided (chair conhguration in Fig. 1(a)) and the single-sided hydrogenations 
(top conhguration in Fig. 1(b)) are taken into consideration to reveal that there are two 
basic types of electronic structures, which could be further modulated by various hydrogen 
concentrations. This work shows many kinds of energy dispersions, including the presence 
or the absence of a Dirac cone, the H-Si-bond-induced partially hat bands, and the nearly 
dispersionless band at the Fermi level induced by broken mirror symmetry. All the pre¬ 
dicted band structures are further rehected in various feature-rich densities of states (DOS), 
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such as strong H-concentration-related peaks. The effect of hydrogenation on band struc¬ 
tures and DOS could be understood by a thorough analysis of the charge distribution of vr 
and a bondings. The predictions of the band gaps and dispersions, and critical character¬ 
istics in DOS could be further identified by ARPES and STS,^’^’^^’^® respectively, with the 
latter being very useful in determining the hydrogen concentration due to its relation with 
the peak intensities. Hopefully, these rich fundamental features in hydrogenated silicenes 
structures can promote potential applications in energy materials or electronic devices. 

2. Methods 

Our first-principle calculations are based on the density functional theory (DFT) imple¬ 
mented by the Vienna ab initio simulation package (VASP).^° The local density approxima¬ 
tion (PDA) is chosen for the DFT calculations, and the projector augmented wave (PAW) 
is utilized to describe the electron-ion interactions. A vacuum space of 10 angstroms is 
inserted between periodic images to avoid interactions, and the cutoff energies of the wave 
function expanded by plane waves were chosen to be 500 eV. For calculations of the elec¬ 
tronic properties and the optimal geometric structures, the first Brillouin zones are sampled 
by 100 X 100 X 1 and 12 x 12 x 1 A:-points via the Monkhorst-Pack scheme. The convergence 
of the Helmann-Feymann force is set to 0.01 eV A ^. 

3. Results and Discussion 

There exist specific relationships between the optimal geometric structures of hydro¬ 
genated silicene and various hydrogen conhgurations. Two different atomic configurations 
of fully hydrogenated silicene are taken into account. The chair-like configuration, where 
H atoms alternatively appear on both sides of silicene plane is shown in Fig. 1(a). The H 
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atoms are set to be in a uniform distribution with each H atom being separated from its 
six nearest H atoms by the same distance (Fig. 1(c)), and the back-site H atoms (green 
circles) are located at the middle of three front-site H atoms (red circles). The hydrogen 
concentration is reduced from 100 % to 25 %, and then 6.3 %, in which the hydrogen- 
silicon ratios correspond to 2 : 2, 2 : 8, and 2 : 32, respectively (Fig. 1(c)). For the top 
configuration, hydrogen atoms only appear on a single side of the plane (Fig. 1(b)), i.e., 
all the back-site H atoms are removed, so the H-concentration is half of the chair-like one. 
These two types are revealed to have similar geometric deformations; however, they may 
lead to diverse electronic properties, mainly owing to the breaking of the mirror symmetry. 

The main characteristics of geometric structures, including the H-Si-Si angles, Si-Si 
bond lengths, and H-Si bond lengths, are dominated by the H-concentration, as shown 
in Table 1. The H-Si-Si angles present a signihcant change from 107.0° to 111.7°, as the 
H-concentration is reduces from 100 % to 3.1 %. The Si-Si bond lengths are non-uniform 
in hydrogenated silicene, while they are at the same value of 2.25 A in pristine silicene. 
The strong H-Si covalent bonds are formed by the obvious attraction of electrons from Si 
atoms to H (discussed later). Such charge transfer makes the nearest Si-Si bond weaker 
but the second-nearest Si-Si bond stronger, as view from the H atom. As a result, the bond 
lengths get longer and shorter, respectively, in the range of 2.31 2.33 A and 2.23 ~ 2.24 

A. It should be noted that the planar Si-Si bond lengths are mainly determined by the a 
bondings due to the (3s, 3pa;, 3pj^) orbitals. However, the H-Si bond lengths are almost 
insensitive to the change in H-concentration, in which they are kept at 1.50 A because of 
the strong and stable covalent a bondings. The above-mentioned bonding configurations 
are closely related to the strong covalent bonds between dp^ and Is orbitals, and the weak 
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sp^ hybridization of four orbitals (3s, 3pj^, 3p2) with 3p2 orbitals perpendicular to the 

plane (discussed later in Fig. 4(a)). These rich geometric configurations lead to the diverse 
electronic properties. 

Hydrogenated silicene exhibits very special energy bands, being enriched by the com¬ 
petition of strong H-Si bondings and weak sp^ hybridization. Pristine silicene, as shown in 
Fig. 2(a), has an almost vanished gap of ~ 1 meV at the K point, and its low-lying linear 
dispersions become parabolic bands with a saddle point at the M point. Such bands can 
be regarded as tt bands, mainly coming from the bonding of 3p^ orbitals of two nearest Si 
atoms. The weak sp^ hybridization indicates a clear separation of the vr bands and a bands 
within ±2 eV. The occupied and unoccupied a bands, with parabolic energy dispersions, 
are initiated at about —1.2 eV and 1 eV, respectively. The highest occupied states (HOS) 
are doubly degenerate at the F point, while the lowest unoccupied states (LUS) are non¬ 
degenerate at the M point. The former energy bands are closely related to the (dp^,, 3py) 
orbitals with their orbital hybridization with dp^ orbitals at lower energy (~ —2.5 eV). The 
third a band with an initial energy at —11.2 eV, dominated by the 3s orbitals (not shown), 
gradually increases with a parabolic behavior. However, it becomes partially flat at —3.2 
eV near the F point, where the 3s orbitals hybridize with the dp^ orbitals (discussed later). 
The main features of silicene are dramatically changed after hydrogenation. 

The chair-like systems exhibit semiconducting behavior with an H-concentration-dependent 
energy gap. Fully hydrogenated silicene does not exhibit a linear Dirac cone as a result of 
the strong H-Si bonding, as illustrated in Fig. 2(b). It has a large indirect band gap of 
Eg = 2.1 eV, which is determined by the HOS and the LUS at the F and the M points, 
respectively. All 3p2 orbitals are saturated by Is orbitals so that the low-lying tt bands 
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are replaced by the a bands corresponding to Si atoms, with the amount of contribution 
indicated by the radius of red circles in Fig. 2(b). Compared to Fig. 2(a), a clear revelation 
is that the valence a bands are formed by the (3pa;, 3py) orbitals. On the other hand, a 
weakly dispersed band at about -3.5 eV mainly originates from one of the particular Is 
orbital of two H atoms (blue circles), and the other Is orbitals contributes equally to two 
partially flat bands. The former characteristic is very useful in gauging the H-concentraion. 
Such electronic states could be roughly regarded as localized states. The H-dependent lo¬ 
calization behavior reflects the fact that the Is orbital has much lower bound state energy 
than the other orbitals. Moreover, the hybridization of Is and dp^ orbitals puts the vr bands 
far away from Ep = 0, and thus the energy gap is determined by the a bands. 

As the concentration decreases, the H-Si bonds no longer dominate, but compete with 
the effect of weak sp^ hybridization to determine the band structures. The 25 % and 6.3 % 
systems, namely the intermediate states between fully hydrogenated and pristine silicene, 
respectively, possess 1.3 and 0.2 eV direct band gaps at the F point, as shown in Fig. 2(c) 
and 2(d). Their two pairs of energy bands nearest to Ep = 0 exhibit weak dispersions with 
narrow band widths, in which they originate from the dp^ orbitals of non-passivated silicon 
atoms. This indicates that the low-lying bands are gradually changed from a bands to vr 

bands. For the middle-energy states (-1 eV-3.5 eV), the a bands mainly come from 

both the passivated (red circles in Fig. 2(c)) and non-passivated silicon atoms. The latter 
contribute much to such bands in the 6.3 % systems as a result of the higher percentage. 
The weakly dispersed band mainly contributed by hydrogen atoms (blue circles) moves to 
-4.0 and -4.2 eV for 25 % and 6.3 % systems, respectively. For other low concentrations, the 
H-dependent energy gaps, as shown in Fig. 2(e), start to fluctuat and gradually decline. 
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As expected, the tunable direct band gap in hydrogenated silicene can be used to provide 
potential applications in nanoelectronic and nanophotonic devices. 

In the top conhguration, the destruction of the mirror symmetry induces strong localized 
states at the Fermi level. The 50 % system has a nearly dispersionless energy band at Ep = 
0 along different directions in Fig. 3(a), which is only associated with three non-passivated 
silicon B atoms (green circles) nearest to the passivated silicon A atom. Such quasi-zero- 
dimensional electronic states are hardly affected by the H-concentration, as indicated by 
the 12.5 % and 3.1 % systems in Figs. 3(b) and 3(c), respectively. The tt bondings between 
B atoms and the A atoms are severely suppressed by the large charge transfer between 
the latter and the H atoms. As a result, B atoms cannot form the extended vr-electronic 
states through interactions with A atoms. It should be noted that these systems are 
gapless, while the conducting behavior is expected to be forbidden by the almost vanishing 
group velocities. The critical difference between the top and the chair conhgurations is 
the intriguing localized states. In addition, the reduction of one effective Is orbital in top 
systems denies the presence of the weakly dispersed band of the chair conhgurations. 

There exist certain similarities between the top and chair conhgurations, including the 
nucleation of tt bands at low H-concentration, the Si-Si bond-related a bands, and the H-Si 
bond-induced localized states. As the H-concentration decreases, both systems exhibit the 
Si-atom-contributed low-lying bands with stronger dispersions and larger energy widths, 
as shown by the 6.3 % (chair) and 3.1 % (top) systems in Figs. 2(d) and 3(c), respectively. 
Such bands mainly result from dp^ orbitals, meaning that the tt bondings are reformed 
among many of the non-passivated Si atoms. For the a bands within -1 —3.5 eV (Figs. 
2(a)-2(c) and 3(a)-3(c)), the initial energy, band width, and state degeneracy are much alike 


between two configurations, owing to the barely infiuenced planar a bonds. Moreover, the 
partially fiat bands related to H-Si bonds (blue and red circles) are confined in a similar 
energy range compared to their counterparts in the chair configuration, refiecting the low 
bound state energy of the H atoms. The above-mentioned similarities are hardly affected 
by the breaking of the mirror symmetry. 

The main characteristics of the energy bands can be verified by ARPES measurements, 
as done for silicenes and few-layer graphenes.Such experiments have been utilized to 
confirm the graphene-like linear energy dispersions for a silicene layer on Ag(lll),^ and 
also to observe the opening of a substantial energy gap in hydrogenated graphene. 

As expected, in the verification of hydrogenated silicene, ARPES can directly detect the 
gradually recovered vr bands, H-related partially fiat bands, and nearly dispersionless bands 
at Ep = 0, which are, respectively, induced by the strengthened tt bonding between Sp^ 
orbitals, the strong H-Si bonds, and the broken mirror symmetry. It is deduced that 
the comparison between theoretical and experimental results of band gaps and energy 
dispersions can provide valuable information in determining the H-concentrations. 

The carrier density (p) and the variation of carrier density (Ap) can provide very useful 
information in the charge distribution in bondings of all orbitals and thus explain the 
dramatic change of low-lying energy bands. The former directly reveals the bond strength 
of Si-Si and H-Si bonds, as illustrated in Figs. 4(a)-4(d). Between two Si atoms of pristine 
silicene or two non-passivated Si atoms of hydrogenated silicene, p shows a strong covalent 
a bond, as indicated by an enclosed black rectangle region in Figs. 4(a), 4(c), and 4(d). 
Such bonds become a bit weaker when the Si atom bonds with an H atom (gray rectangle 
in Figs. 4(b)-4(d)). The little change is responsible for the similarities of a bands between 
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the two configurations under various H-concentrations. Away from the central line between 
two Si atoms, carrier density is lower, but its distribution is extended between them (purple 
rectangle). This is evidence for the vr-electronic states, revealing the gradual recovery of 
TT bands at a low H-concentration. However, barely extended states, corresponding to the 
TT-electron-depleted region, are observed between the passivated and non-passivated atoms 
(red rectangle in Fig. 4(d)). Such charge distribution is similar to that of isolated Si 
atoms and demonstrates the disappearance of vr bonding. This also induces the nearly 
dispersionless band in all top configurations. On the other hand, the charge density of 
H-Si bonds (dashed rectangle) is much stronger than those of all Si-Si bonds. The very 
strong bond strengths are insensitive to H-concentration in both two systems, so that the 
H-related energy bands are confined at about -3.5 eV —5.5 eV. It is also noted that the 
charge distributions can explain why only the nearest Si-Si bond lengths are changed, as 
discussed earlier in Table I. 

In order to comprehend the charge transfers among all orbitals, the detailed variation 
of carrier density is clearly illustrated in Figs. 4(e)-4(h). Ap is created by subtracting the 
carrier density of isolated silicon atoms (isolated silicon and hydrogen atoms) from that of 
a silicene (hydrogenated silicene). For pristine silicene, as indicated in Fig. 4(e), electrons 
are transferred from the blue region near Si atoms to the location at the middle of two 
Si atoms (yellow region within a black rectangle), forming covalent a bonds. Distantly 
away from the middle region between two Si atoms, p also increases in the green region 
(purple rectangle), which verifies the creation of vr bondings. On the other hand, all the 
hydrogenated cases display high charge transfer from Si atoms to H atoms (~ 0.9 e by using 
Badar charge analysis), forming strong H-Si bonds (dashed rectangle Figs. 4(f)-4(h)). This 
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also accounts for a rather weak bonding strength of the nearest Si-Si bond around the H 


atom. 

The density of states (DOS) directly reflects the primary characteristics of the band 
structures, in which the partial DOS (PDOS) can be used to fully understand the orbital- 
decomposed contributions. For pristine silicene, DOS is zero at Ep = 0; it possesses a 
linear E-dependence at low energy, and then a symmetric logarithmic divergent peak at 
E ~ ±1 eV (black curve in Fig. 5(a)). Such features, respectively originating from the 
almost vanished energy gaps, isotropic linear bands, and the saddle points of vr bands, are 
dominated by the Sp^ orbitals (purple curve). The middle-energy DOS (-1 eV —3 eV) 
reveals a shoulder (-1.2 eV) structure and a prominent peak (-2.7 eV), corresponding to 
the maximum point and the saddle M point of the a bands, respectively. They are mainly 
contributed by the Sp^, and 3py orbitals (green and blue curves). Also, the hybridization 
between dp^ and 3s orbitals induces an extra strong peak at -3.2 eV. These characteristics 
further illustrate that vr and a bands have the weak interactions at low energy. 

The low energy DOS is drastically changed after hydrogenation. For the fully hydro¬ 
genated silicene in the chair conhguration, DOS displays a vacant region within E ~ ±1 
eV, and an extra strong peak at about -3.5 eV and -4.5 eV, that comes from the H-related 
weakly dispersed band and the partially flat band, respectively, as shown in Fig. 5(b). 
These features are mainly owed to the strong bonding of dp^ and Is orbitals (purple and 
red curves). However, the shoulder structure at —1.0 eV and the prominent peak at —2.5 
eV, constituted by the 3pa; and 3pj/ orbitals, are very similar to that of pristine silicene as a 
result of the small change in the a band. As the H-concentration decreases, two prominent 
peaks are revealed at ±0.6 eV in the 25 % system, as shown in Fig. 5(c). These peaks 
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have widths of about 0.5 eV, and their states are associated with vr bands that are related 
to Sp^ orbitals. The middle energy states are also dominated by the 3pa; and 3pj, orbitals; 
however, the single peak structure in the 100 % system changes into several subpeaks due 
to an increase of the number of subbands. In the 6.3 % system, 3p2 orbitals contribute most 
of the electronic states within ±1 eV, indicating a similar tt band width to pristine silicene. 
The two lower concentrations show two strong peaks between -4.0 —4.5 eV, originating 
mainly from the Is orbital of the H atoms (red curve) and the Op^ orbitals of Si atoms 
(purple curve). Nevertheless, the intensities of these two peaks are roughly proportional to 
the H-concentration, especially for the peak closer to the Fermi level. 

Common features in DOS between the top and chair systems are the strong peaks 
associated with the planar covalent a bonds and the H-related bonds. Analogies are 
found between a top case and its corresponding chair conhguration (with double the H- 
concentration), including the widened low-energy vr band width (Figs. 5(d) and 5(g)), 
the two 3p2-orbital-induced prominent peaks at ±0.6 eV (Figs. 5(c) and 5(f)), the dp^, 
and 3py orbitals dominating middle-energy peaks (Figs. 5(b)-5(g)), and the strong peaks 
contributed by the H-Si bonds. Such special peaks are located in almost the same energy 
range, while their intensities are gradually decreased with diminishing H-concentration. 
On the other hand, the destruction of the mirror symmetry leads to a critical difference, 
namely a delta-funtion-like peak at Ep (Figs. 5(e)-(g)). Its intensity gradually decreases, 
directly reflecting the percentage of the nearest non-passivated Si atoms. Both dp^ and 
3s orbitals contribute to these peaks, as shown by the purple and red curves, respectively. 
This directly illustrates that the non-extended vr-electrons of dp^ orbitals are hybridized 
with 3s orbitals. 
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The STS measurements, in which the tunneling differential conductance map of the 
dl/dV-V curve is proportional to DOS, can provide an accurate and efficient way to examine 
the H-concentration. It has been adopted to identify the Dirac point of pristine siliceneT 
Hydrogenated graphene that possesses a gap larger than 2 eV has also been measured with 
this methodff® The STS measurements could be used to determine the energy gaps and the 
specific DOS peaks of hydrogenated silicene, especially the strong H-bond-related peaks 
(red triangles in Figs. 5(b)-5(d)), and the delta-funtion-like peak at Ep (Figs. 5(e)-5(g)). 
These intensity of these peaks is predicted to be strongly associated with the percentage 
of H atoms. The intensity of the H-related peaks can be well fit by a linear relationship 
I = 1.6Io{N/No), where I is the peak intensity and N is the H-concentration, with Jq and 
Nq being the value for the highest concentration. For the delta-funtion-like peak at Ep 
this relation changes into I = 1.3Iq{N/Nq). Given the linear dependency of DOS peaks 
associated with the contributions from H atoms and the nearest non-passivated atoms, we 
have a feasible way to identify the hydrogen concentration. 

4. Conclusion 

The geometric and electronic properties of hydrogenated silicene are investigated by 
DFT calculations. Chemical and physical properties are enriched by different hydrogen 
configurations and concentrations, in which the special H-Si bonds are the critical factor 
affecting the optimized geometric structures, charge distributions, energy bands, and DOS. 
The high charge transfer between Sp^ and Is orbitals, directly or indirectly leading to 
variations in the bond lengths and carrier densities, is responsible for the observed dramatic 
change in the electronic properties. The three major revelations of the band structures are 
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highlighted by the absence or recovery of low-lying vr bands, the H-related partially flat 
bands, and the nearly dispersionless band at Ep = 0. The first one is determined by the 
competition of strong H-Si bondings and weak sp^ hybridization. The band gap mainly 
determined by the a bands is largest for the 100 % H-concentration. As the concentration 
drops, the influence of the vr bands gradually rises, causing the gap size to shrink. The 
second one, purely originating from the strong charge bondings between H and Si atoms, 
is a common characteristic for both top and chair systems under various H-concentrations, 
while the last one, being ascribed to the destruction of mirror symmetry, is a distinguished 
characteristic for top configurations. As for DOS, the recovery of vr-bands-induced low- 
lying peaks, the diminishing of H-related strong peaks, and the weakened delta-funtion- 
like peak at Ep directly reflect the reduced H-concentration. The feature-rich band gaps 
and energy dispersions, and the strong DOS peaks, respectively, could be measured by 
ARPES and STS experiments, giving useful information about the hydrogen distribution 
and concentration. Especially, the linear dependence of the special DOS peaks provides 
an effective path to accurately detect the H-concentration. Hydrogenated silicene is very 
suitable for the application of hydrogen storage, and the tunable electronic properties might 
be potentially important for applications in nanoelectronic and nanophotonic devices. 
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Figure 1: Side view atomic configurations of hydrogenated silicenes: (a) chair conhguration, 
where the H atoms appear on both sides of silicene plane, and (b) top conhguration, 
where all H atoms are located on the same side of the plane, (c) top view for various 
concentrations, namely 2:2, 2:8, and 2:32. 21 
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Figure 2: Energy bands of (a) pristine silicence, and hydrogenated silicenes with chair 
configuration for (b) 2:2, (c) 2:8, and (d) 2:32 hydrogen coverages. The inset of (a) shows 
the low-lying energy bands, and (b) the band gaps of various concentrations. 
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Figure 3: Energy bands of top configurations for (a) 1:2, (b) 1:8, and (c) 1:32 concentrations. 


























































Figure 4: The carrier density p of (a) pristine silicene, and hydrogenated silicenes with (b) 
2:2, (c) 1:2, and (d) 1:8 concentrations. The variation of carrier density Ap of (e) pristine 
silicene, and hydrogenated silicenes with (f) 2:2, (g) 1:2, and (h) 1:8 concentrations. 
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Figure 5: Densities of states of (a) pristine silicenes, chair configurations with (b) 2:2, (c) 
2:8, and (d) 2:32 hydrogen coverage, and top conhgurations with (e) 1:2, (f) 1:8, and (g) 
1:32 hydrogen coverage. 
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